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Association of genital human papillomavirus 
infection with HIV acquisition: a systematic review 
and meta-analysis 

Pascale Lissouba/ Philippe Van de Perre,^'^'^ Bertran Auvert^'^'^'^ 



ABSTRACT 

Objectives To conduct a systematic review and a 
meta-analysis of epidemiological studies investigating the 
association of genital human papillomavirus (HPV) 
infection and HIV acquisition. 
Design Systematic review and meta-analysis. 
Data Sources Scientific databases and conference 
abstracts were systematically searched to identify all 
relevant studies published up to 31 January 2012. 
Search terms included 'HIV, 'HPV, 'human 
papillomavirus' and 'papillomaviridae' as keywords or 
text, in the title or abstract. 
Methods To be eligible for inclusion, a study had to 
be conducted among humans, report data on HIV 
incidence, and assess genital HPV infection. Summary 
ORs and 95% CIs were estimated from the extracted 
data using random-effect meta-analysis. Subgroup 
analyses were conducted for high-risk (HR) and low-risk 
(LR) HPV oncogenic risk groups. Between-study 
heterogeneity and publication bias were assessed. 
Results Of 2601 identified abstracts, six observational 
studies, comprising 6567 participants were retained for 
the systematic review and the meta-analysis. HIV 
acquisition was significantly associated with HPV 
infection (summary OR=1.96; 95% CI 1.55 to 2.49). 
HIV incident infection was significantly associated with 
HR-HPV in five of six studies and with LR-HPV in two 
out of five. The association was significant for HR-HPV 
(summary OR=1.92; 95% CI 1.49 to 2.46) and 
borderline for LR-HPV. No between-study heterogeneity 
was detected. There was a borderline indication of 
publication bias. 

Conclusions Further research is needed to elucidate 
the biological mechanisms involved, and assess the 
effect of HPV vaccination on HIV acquisition, using 
vaccines with broad coverage of HPV genotypes. Such 
research could have important public health implications 
for HIV prevention. 



INTRODUCTION 

Evidence from systematic reviews and meta-analyses 
of epidemiological studies suggests that sexually 
transmitted infections (STI) such as Neisseria gonor- 
rhoea, Chlamydia trachomatis, Trichomonas vagina- 
lis, infectious syphilis and herpes simplex virus-2 
(HSV-2) play a role in HIV acquisition. ^""^ A model- 
ling study has estimated that over 50% of HIV 
sexual transmission could be attributed to STI.^ The 
contribution of STI to the HIV epidemic is still 
being debated, due to the observational nature of 
the evidence,^ and the poor efficacy of randomised 
clinical trials of STI management on HIV 



acquisition/ Nonetheless, there is biological plausi- 
bility to the causal association of STI Wixh HIV 
acquisition, w^hich is corroborated by the enhanced 
susceptibility of STI-affected individuals to HIV^ ^ 

Little research had been conducted on the associ- 
ation of HIV acquisition With, genital human papillo- 
mavirus (HPV) infection, the most common STI 
M^orldw^ide,^^ and a major cause of morbidity and 
mortality, contributing to up to 5.2% of the cancer 
burden in the w^orld.^^ Genital HPX w^hich targets 
the basal epithelial cells of the mucosal surface,^^ is 
highly infectious, w^ith transmissibility estimated to 
be higher than any other viral STI.^^ Mucosal HPV 
genotypes are divided into high-risk (HR) and 
loM^-risk (LR) HPV based on their association w^ith 
cervical malignancy. HR-HPV genotypes cause cer- 
vical carcinoma, one of the leading causes of cancer 
deaths in the w^orld, and have been implicated in 
vulvar, vaginal, penile and anal cancers. LR-HPV 
genotypes are associated With, genital Mparts, w^hich 
negatively impact quality of life and psychosocial 
functioning.^^ HPV prevalence is estimated to 
range from 1.6% to 41.9% among w^omen,^^ and 
1.0% to 95.0% among men,^'^ according to factors 
such as setting, age and sexual behaviour. 

To date, the bulk of research has investigated the 
link between HPV infection and HIV-serostatus. 
Studies indicate that HIV seropositivity is strongly 
associated w^ith higher HPV prevalence, higher 
HPV incidence,^^ a reduced likelihood of HPV 
infection clearance w^hich increases the risk of 
related lesions,^^ and higher prevalence of 

infection w^ith multiple HPV genotypes.^^ 
HIV seropositivity has also been linked to high- 
grade and low^-grade cervical squamous intraepithe- 
lial lesions,^^ as w^ell as invasive cervical 
carcinoma, a recognised AIDS-defining condition^^ 
and a common cancer among HIV-positive w^omen 
in sub-Saharan Africa.^^ 

Research on the role of genital HPV in HIV inci- 
dent infection is currently limited, w^hereas it is bio- 
logically plausible. HPV could predispose to HIV 
infection and dissemination through its disruption 
of the epithelium integrity and the mucosal 
immune system,^'^ by (A) altering the density and 
functional activities of epithelial Langerhans cells 
w^hich capture infecting pathogens,^ (B) enabling 
the recruitment and activation of HIV target cells 
such as T-lymphocytes,^^ (C) reducing the expres- 
sion of proteins acting in cell adhesion and tumour 
suppression,^^ (D) dow^nregulating proteins that 
promote the infiltration of Langerhans cells 
through the epithelium,^ ^ (E) reducing the 
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production of proteins involved in antimicrobial activities'^ and 
(F) upregulating inflammatory cytokines,'^ w^hich may 
increase HIV replication. 

There have been several calls for the investigation of the associ- 
ation of genital HPV infection w^ith HIV acquisition, dating as far 
back as 1993/ In view of the scope and severity of 

both infections and their synergy, a review^ of the existing epi- 
demiological evidence is w^arranted. Since there are now^ two 
highly efficacious vaccines against some HPV genotypes,^^ this 
research may have important implications for policies on HPV 
vaccination and public health efforts to limit the spread of HIV 

OBJECTIVES 

The objectives of this study w^ere to systematically review^, sum- 
marise, and conduct a meta-analysis of epidemiological studies 
investigating the association of genital HPV infection w^ith HIV 
incident infection. 

METHODS 

The review^ and meta-analyses w^ere performed according to the 
PRISMA 2009 (Preferred Reporting Items for Systematic 
Review^s and Meta-Analyses) and STROBE (Strengthening the 
Reporting of Observational Studies in Epidemiology) 
guidelines.^^ 

Search strategy and selection criteria 

A bibliographic search w^as performed in MEDLINE, EMBASE, 
Web of Science and the Cochrane database, to identify all rele- 
vant epidemiologic studies pubHshed up to 31 January 2012. 
Search terms consisted of ('HIV OR 'human immunodeficiency 
virus') AND 'HPV OR 'human papillomavirus' OR 'papilloma- 
viridae') as keyw^ords or text, in the title or abstract. There w^ere 
no language restrictions. This broad search strategy ensured that 
all preidentified eligible studies were included in the results. 
Abstracts from relevant conferences w^ere also searched. 

All retrieved abstracts w^ere read by two of the review^ authors 
(PL and BA), and evaluated for eligibility. To coincide w^ith the 
study objectives and be selected, a study had to be conducted 
among humans, report data on HIV incidence, and assess its 
association w^ith genital HPV infection. Genital HPV could be 
reported either by oncogenic risk group (LR-HPV vs HR-HPV), 
individual genotypes, or as a categorical (prevalent HPV infec- 
tion), or discrete (number of HPV genotypes) variable. To assess 
HPV as a potential risk factor of HIV acquisition, only studies 
assessing HPV infection before measuring HIV incidence w^ere 
retained. Cohorts could include HIV seropositive individuals at 
baseline. Review^s, cross-sectional studies, case studies and case 
series w^ere excluded. 

Data extraction 

The full text of the retained articles wsls independently read by 
all review^ authors. Data were extracted using a standardised 
form. Information obtained included HPV prevalence in the 
sample, number of HIV incident infections, unadjusted and 
adjusted estimates, including 95% CIs, of the association 
between HPV infection and HIV acquisition, type of estimate, 
study setting, sample size, demographic characteristics, HPV 
sampling methods, length of follow^-up and covariates used for 
adjustment in the multivariate models. 

Statistical analyses 

The summary ORs and 95% CI w^ere estimated using 
random-effect meta-analytic techniques because there w^as sub- 
stantial heterogeneity betw^een studies in terms of sample size. 



participants characteristics, reported HIV incidence, length of 
follow^-up, covariates used for adjustment and HPV assessment 
(including sampling method, oncogenic risk group classification 
and measured prevalence). 

The main analyses were conducted for HPV infection, regard- 
less of oncogenic risk group. Adjusted estimates were used 
M^henever possible. Sub-group analyses w^ere conducted by HPV 
oncogenic risk group. One study^^ reported estimates for three 
different timings of HPV assessment: at baseline, concurrently 
w^ith HIV testing and w^ithin 6 months of HIV testing. This 
latter assessment, labelled 'recent' HPX was used in the main 
analysis. Sensitivity analyses were conducted w^ith the other two 
timings of assessment. 

For each summary analysis, betw^een-study heterogeneity w^as 
assessed using the I-squared statistic,^^ w^hich quantifies the 
degree of heterogeneity. Publication bias was evaluated using the 
Begg-adjusted rank correlation test, w^hich assesses the associ- 
ation between the treatment effect and its variance,^^ and the 
Egger regression asymmetry test.^^ All statistical analyses w^ere 
performed using the Stata Statistical Software: Release VI 0 
(StataCorp 2007, College Station, Texas, USA: StataCorp LP). 

RESULTS 
Study selection 

The database search identified a total of 2601 abstracts (2170 
articles and 431 conference proceedings). The study selection 
process is detailed in figure 1. After the exclusion of duplicate 
citations and non-relevant abstracts, 14 published studies w^ere 
retained. Eight studies were rejected because they did not fulfil 
the inclusion criteria. Live nested cohort studies, and one case- 
control study on the association betw^een genital HPV infection 
and HIV acquisition w^ere identified."^^ 42^6 ^eie^;^^^^ studies are 
summarised in table 1. 

Study characteristics 

Eour studies were conducted among w^omen^^ and two 

among men."^^ They included a total of 6567 participants. All 
but one study among adult men who have sex w^ith men^^ w^ere 
conducted in sub-Saharan Africa among heterosexual adults. 
Tw^o of the African studies w^ere conducted among high-risk 
M^omen who were mostly sex w^orkers.^^ 

HPV prevalence across studies reporting baseline overall HPV 
infection ranged from 24.5% to 56.8%. In studies conducted 
among w^omen, HPV testing w^as performed on cervical 
samples. HPV sampling among men included the collection of 
anal samples in the Chin-Hong and colleagues study^^ and 
exfoliated cells from the glans or coronal sulcus and the penile 
shaft in the Smith and colleagues study.^^ HPV sampling was 
performed by health professionals, apart from the Smith- 
McCune and colleagues study^^ w^hich included self-collected 
and clinician-collected cervical sw^abs. 

In all six studies, HPV testing w^as conducted using PCR and 
positive HPV results were genotyped using type-specific probes 
on PCR products. How^ever, these methods differed in terms of 
HPV genotypes detected and their classification by oncogenic 
risk group, w^ith some genotypes classified as high-risk in some 
studies and as low^-risk in others (see online supplementary 
appendix table S2a and S2b). Six studies reported estimates for 
HR-HPV and five for LR-HPV In the subgroup analyses, the dis- 
tinction between HR-HPV and LR-HPV w^as made according to 
the criteria reported in each study because genotype-specific 
data w^ere not available. 

Time betw^een HPV status and HIV seroconversion assess- 
ments varied betw^een studies. When HPV was assessed at the 
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Figure 1 Flow diagram of study 
selection for the review and 
meta-analysis. Of the 2601 abstracts 
retrieved, 14 were selected for full 
reading. Based on the selection 
criteria, six studies were retained for 
the systematic review and the 
meta-analyses. 



Abstracts identified through database 
searching 
(2170 published + 43 1 abstracts) 
n = 2601 



Abstracts screened 
n = 2601 



Abstracts excluded for non- 
relevance 
n =2587 



Full-text articles assessed for 
eligibility 
n=14 



Studies included in 
qualitative synthesis and 
meta-analysis 
n = 6 



Full-text articles excluded 
(n = 8) 

Reasons : 

Improper HPV assessment 

(n-6) 

No HIV incidence data (n=l) 
Insufficient data to derive an 
estimate (n=l) 



start of the follow-up period over which HIV incidence was 
measured (or at month 6 for one study"^^), follow-up varied 
between 21 and 42 months.^^ The Smith-McCune 

and colleagues study also assessed incident HIV infection 
detected concurrently and within 6 months of HPV infection 
assessment. The Averbach and colleagues study^"^ used HPV data 
detected on average 80 days before HIV seroconversion. 

HPV genotypes found to be associated with HIV acquisition 
included HPV 52,^^ as well as HPV 31, 58 and 70.^^ In terms 
of the four HPV genotypes covered by the currently available 
vaccines, only one study investigating individual HPV genotypes 
found significant positive associations between HIV acquisition 
and genotypes 16 and 18.^^ None of the studies found associa- 
tions with genotypes 6 and 11. 

The Auvert and colleagues study,"^^ the Averbach and collea- 
gues study,^^ and the Chin-Hong and colleagues study,^^ found 
that the risk of HIV acquisition increased with the number of 
detected HPV genotypes. By contrast, in the Smith and collea- 
gues study,^^ risk of incident HIV infection did not differ 
between individuals with multiple infections versus those 
infected with only one genotype. 

Both the Smith-McCune and colleagues study"^^ and the 
Averbach and colleagues studies^^ investigated HPV persistence 
and clearance, and found a significant association between non- 
persistent HPV and HIV acquisition. However, the definition of 
non-persistence differed between the studies. 



Study quality 

For all six studies, unadjusted and adjusted estimates, as applic- 
able, and their 95% CI are reported in table 1. All but one 
study^^ reported adjusted OR. Statistically significant estimates 
are reported in bold. Variables included in the multivariate 
models are listed in onHne supplementary appendix table S3. 
Most studies controlled for age "^"^-^^ and HSV-2 

status.^ ^ All studies conducted among women, and 

reporting adjusted estimates, controlled for condom use and 
sexual behaviour.^ ^ Two studies conducted in sub-Saharan 
Africa^ ^ controlled for circumcision status, known to be asso- 
ciated with HIV incidence. Measures of sexual partners' charac- 
teristics and participants' sexual behaviour were 
population-specific, assessed at different times during follow-up, 
and varied substantially across studies. 

Main meta-analysis 

Six studies were included in the meta-analysis.^^ 42^6 
reported significant estimates of the association of HPV with 
HIV incidence.^^ "^^^^ As shown in figure 2, the computed 
summary OR indicated that overall, individuals with genital 
HPV infection, regardless of oncogenic risk group, had twice 
the risk of acquiring HIV (summary OR, 1.96; 95% CI, 1.55 to 
2.49). There was no evidence of between-study heterogeneity 
(p = 0.484), suggesting that a pooled analysis was suitable. 
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Table 1 Characteristics of identified studies investigating the association between genital HPV infection and HIV acquisition 
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3.6 (1.5 to 8.4)* 
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2009 




Men who have 
sex with men 


older 












HR-HPV: 26.4% 
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NR 
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Smith-McCune 
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1987 
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21 


88 


HR 


Baseline HPV: 


NR 
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et al 2010 
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12-month visit; quarterly vaginal 
swabs 








24.5% 
Baseline 
HR-HPV: 16.1% 
Baseline 
LR-HPV: 10.8% 


NR 
NR 


2.0 (1.2 to 3.2) 

1.5 (0.92 to 2.4) 


Smith etal 2010 


Kenya 


Men general 
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18-24 
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Exfoliated penile cells from glans/ 
coronal sulcust at baseline 


42 


63 


HR 
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HR-HPV*: 
34.8% 
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15.5% 
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NR 
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1.9 (0.9 to 3.6) 


Veldhuijzen et al 
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Cervical samples at month 6 
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OR 


HPV: 47.8% 


NR 


NR 
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HR-HPV: 33.3% 
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Women 


median 


88 


first available cervico-vaginal 
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25 


HR 


HPV: 


NR 
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High-risk (all sex 
workers) 


24 




rinses of 3 obtained during 
follow-up, before seroconversion 








HR-HPV: 70.5% 
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1.7 (1.0 to 2.7)§ 
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Case-control Study 
























Averbach ef al 


Zimbabwe 


Women general 


18-35 


1 45 cases 


Cervical swabs at t_i (to=time of 


to-t_i=80.5 days 
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OR 


HPV: 49.1% 


2.7 (1.7 to 4.1) 
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HIV detection) 
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Any HR-HPV: 
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12.0% 


2.7 (1.7 to 4.3) 
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2.3 (1.4 to 3.9) 
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Bold indicates statistically significant estimates. 
* Infection with two or more HPV genotypes. 
tSampling penile shaft for HR-HPV and LR-HPV 

*HPV infections with multiple HPV types were considered high risk if one or more HR-HPV type was detected. All the others were considered LR-HPV infections. 
§lnfection with one HPV genotype versus no HPV. 

HPV, human papillomavirus; HR-HPV, high-risk HPV; IRR, incidence rate ratio; LR-HPV, low-risk HPV; NR, not reported; RCT, randomised control trial. 
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Figure 2 Combined estimate of the 
association of human papillomavirus 
infection with HIV acquisition. The 
weights and combined estimate were 
computed using random-effect 
meta-analysis. Estimates from 
individual studies and the combined 
estimate represented on the x-axis are 
in the natural log scale. 
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Sub-group analyses 

When the association with HIV acquisition was investigated by 
subgroups, it was found to be statistically significant with 
HR-HPV (summary OR, 1.92; 95% CI, 1.49 to 2.46; test of 
heterogeneity p = 0.511), and borderline with LR-HPV 
(summary OR, 1.53; 95% CI, 0.96 to 2.42; test of heterogen- 
eity p = 0.029). 

Sensitivity analyses 

The summary ORs computed with the alternative assessments^^ 
for both main and subgroup analyses were almost identical to 
the original estimates (see online supplementary appendix). 

Publication bias 

The funnel plot for the detection of publication bias associated 
with the random-effect meta-analysis of HPV and HIV acquisi- 
tion is presented in figure 3. It appears that studies with larger 
SEs tended to report larger estimates. However, the Begg's test 
was not statistically significant (p = 0.339). Visually, the funnel 
plots seemed asymmetrical, indicating a possible publication bias 
towards studies reporting significant associations. This was con- 
firmed statistically by the Egger's test, which indicated a border- 
line statistically significant publication bias (p = 0.054). 

In the subgroup analyses, no significant publication bias was 
detected for HR-HPV (Begg's test p = 0.806; Egger's test 
p = 0.181). There was a significant publication bias for LR-HPV 
(Begg's test p = 0.089; Egger's test p = 0.022). 

DISCUSSION 

This systematic review and meta-analysis investigates both quali- 
tatively and quantitatively empirical evidence of the association 
between genital HPV infection and HIV acquisition. The fact 
that all studies included have been published recently indicates 
that this is a topic which has been overlooked in past HIV pre- 
vention research. Our findings show that genital HPV infection 
is associated with incident HIV infection, and that this associ- 
ation is significant with HR-HPy and borderline with LR-HPV 



This study is based on evidence from observational research, 
hence it cannot distinguish between causality and confounding. 



Biological mechanisms 

The enabhng effect of HPV on HIV acquisition is plausible, pri- 
marily because HPV is assumed to be more infectious than 
HIX^^ and HPV prevalence exceeds that of HIX even in areas 
where HIV is endemic. Several hypotheses may explain a pos- 
sible higher susceptibility of HPV-infected individuals to Hiy 
including HPV-induced epithelial breaches or displasia, 
T-lymphocytes influx, the production of proinflammatory cyto- 
kines, or a reduced secretion of proteins involved in the 

Begg's funnel plot with pseudo 95% confidence limits 



1 

0.4 

s.e. of: log[or_allr] 



— I — 
0.6 



Figure 3 Funnel plot for the meta-analysis of human papillomavirus 
(HPV) and HIV acquisition The horizontal line represents the natural 
log (In) of the combined OR. The funnel lines represent the 95% 
pseudo confidence limits for this estimate. The graph suggests that 
estimates of the association between HPV and HIV acquisition reported 
by studies with larger SEs tend to be larger. Furthermore, the graph is 
slightly asymmetrical, indicating a possible publication bias towards 
studies reporting significant associations. 
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immunologic response, which may all promote HIV entry and 
replication.^ 

However, the observed association could be due to host suscep- 
tibility for HPV and HIV^^ The two viruses share common mech- 
anical access to susceptible cells and common portals of entry, 
such as the squamo-columnar junctions in the anogenital mucosae, 
which are prime sites of metaplastic conversion, and prone to viral 
establishment.^^ Adhesion molecules are less densely expressed at 
these sites, favouring paracellular permeation.^^ 

Breaches of the mucosal surface, which enable HPV to access 
lower layers of stratified epithelia during its initial invasion 
phase, may favour HIV acquisition as well.^^ 

Lastly, the association could be mediated by other STL For 
instance, HPV is more frequent in HSV-2 infected individuals. 
HSV-2 replication, associated with primary infection and in-situ 
reactivation of latent infections, facilitates HIV acquisition.^^ 
Thus, remote confounding could partly explain the association 
observed. 

Study limitations 

There are several limitations to this meta-analysis. First, it was 
conducted among a small number of studies, so it is underpow- 
ered, which may have affected the I-squared statistic.^^ Second, 
all studies were observational, which increases the likelihood of 
bias. Third, the qualitative review of the studies indicated sub- 
stantial variations in terms of study design, participant character- 
istics, HPV status assessment, and covariates used in the 
adjusted models. Two salient variations were: (A) HPV sampling 
among men, which was not standardised, and could explain the 
differences in HPV and HR-HPV prevalences observed"^^ and 
(B) the classification of HPV genotypes by oncogenic risk group, 
which differed across studies. Since the genotype-specific preva- 
lences were not consistently reported, it was not possible to 
evaluate how these differential classifications may have affected 
the subgroup estimates. Similarly, it was not possible to explore 
genotype-specific associations of HPV with HIV as data was not 
available. Another important issue was the time between HPV 
status assessment and HIV testing, which varied across studies. 
It cannot be excluded that the association is time-dependent, 
and that the transient nature of HPV infection^ ^ could have 
affected its strength. Fourth, metaregressions to examine the 
impact of factors, such as study site, sexual risk behaviour, as 
well as relevant analyses stratified by sampling sites and 
methods, could not be performed for lack of power. Fifth, the 
adjustments for the multivariate estimates used in the 
meta-analyses were performed on different variables. Sixth, pub- 
lication bias was suggested, which could affect the validity of 
the findings. Finally, the generalisation of the study findings is 
difficult since all but one study were conducted in sub-Saharan 
African countries with high incidence of HIV 

Implications 

This study has important implications for future HIV prevention 
research. As a potential risk factor, the contribution of HPV 
infection to the HIV epidemic must be explored comprehen- 
sively, using longitudinal study designs, standardised HPV sam- 
pling and HPV genotypes classification, and analytic methods 
with adjustment on relevant covariates. Further research is 
needed to understand the aetiology of the association, including 
the biological mechanisms that may underlie oncogenic risk 
group-specific or genotype-specific associations. The effect of 
HPV vaccination on HIV acquisition may be assessed since there 
exists two highly effective HPV vaccines which have been 
approved for use among women,^^ and one among men.^^ To 



this day, randomised controlled trials of STI management have 
yielded contradictory results, which may have been partly due 
to the insufficient potency of the intervention, such as in the 
case of HSV-2. ^ As indicated in the present study, there may be 
discrepancies regarding the HPV genotypes for which the vac- 
cines are providing protection against, and those associated with 
HIV seroconversion. Hence, the development and use of vac- 
cines covering a wider range of genotypes is necessary. 
Modelling studies on the potential impact of HPV vaccination 
campaigns on the HIV epidemic are also needed. Such research 
could be an important contribution to the global fight against 
HIV 



CONCLUSION 

The observed association of HPV infection with HIV acquisition 
has important public health implications. If the causal associ- 
ation between HPV and HIV acquisition is verified, HPV vaccin- 
ation, along with other biomedical interventions, has the 
potential to be an additional HIV control strategy, with the 
unique advantage of providing long-term benefits to both men 
and women. 



Key messages 



► The association between HIV acquisition and genital human 
papillomavirus (HPV), a major burden of morbidity and 
mortality worldwide, has received limited attention. 

► It is, however, biologically plausible that HPV enables HIV 
infection and further research is required to elucidate the 
biological mechanisms involved. 

► This review and meta-analysis suggests that HPV-infected 
individuals have twice the risk of acquiring HIV. 

► HPV infection may be a risk factor of HIV acquisition. The 
effect of HPV vaccination on HIV infection should be 
assessed. 
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